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Abstract. The Blandford-Znajek mechanism has long been regarded as a key 
ingredient in models attempting to explain powerful jets in AGNs, quasars, 
blazzars etc. In such mechanism, energy is extracted from a rotating black 
hole and dissipated at a load at far distances. In the current work we 
examine the behaviour of the BZ mechanism with respect to different boundary 
conditions, revealing the mechanism robustness upon variation of these conditions. 
Consequently, this work closes a gap in our understanding of this important 
scenario. 
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1. Introduction 

The Blandford-Znajek (BZ) effect [1] has been proposed as one key mechanism to 
explain energetic jets from black hole systems. In this mechanism a spinning hole 
interacts with magnetic fields sourced by an accretion disk (see. e.g. [31[3J|3]). The 
magnetic field extracts rotational energy from the spinning black hole and dissipates 
this energy at a far away load. The details of this process are complex and its 
understanding requires analyzing the behaviour of the plasma interacting with the 
strong curvature region around the black holes and how the energy extracted, which 
is represented by a powerful collimated Poynting flux, is dissipated at large distances. 

A simple picture explaining this process is provided by the membrane 
paradigm [5], where the system is modeled as a circuit composed of a battery - 
provided by the charge separation induced on the black hole-, two long wires -along 
magnetic field lines- and a load -at large distances from the black hole-. The energy 
tapped from the spinning black hole is understood in terms of an induced EMF 
flowing from the pole to the equator and released at the load. While this picture 
provides a basic understanding of the underlying phenomena, the role of the load in 
the resulting jet and energy dissipation has not been the subject of a detailed analysis. 
In particular, the possibility that different loads could significantly affect the jet has 
not been explored. If such were the case, this mechanism could be questioned as a 
robust ingredient to explain black hole systems producing collimated energy outputs. 

In this note, in order to investigate this issue, we examine the influence of the 
boundary conditions on the dynamics of the system and, in particular, on the resulting 
Poynting flux energy. By considering different conditions, playing the role of the 
load, we illustrate the robustness of this mechanism. Our results indicate that while 
these conditions do result in some differences, these do not significantly affect the jet 
structure and energetics of the system. 

2. Physical system 

The system is modeled by assuming the force free approximation in a curved 
background provided by a Kerr black hole spacetime. The black hole is initially 
embedded in a pure magnetic field configuration, with its intrinsic angular momentum 
(or spin) aligned with the asymptotic magnetic field. We study the dynamics of the 
electromagnetic fields numerically and examine the induced Poynting flux. The details 
of our formulation and implementation are provided in [6j E] ■ We here consider more 
general boundary conditions, which will allow us to examine scenarios corresponding 
to different resistive loads. To define such conditions we must first understand the 
characteristic structure of the system 

2.1. Characteristic decomposition 

Our starting point are the (general relativistic version of the) Maxwell equations. 
Further, we consider the augmented version of the equations, which include the 
divergence-cleaning fields {0, ^} to dynamically control the constraint violations 
[H [10] . This system is defined by 



(d t - Cp)E l - e ijk Vj{aB k ) + a^V^ = atrKE 1 - a.P , 



(1) 
(2) 
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(S t - C p )B i + e ijk V 3 (aE k ) + a^V ' j( j> = atrKB i , (3) 

(S t - Cfs)(j) + aViB 1 = -an<f> . (4) 

It is illustrative to consider first the electrovacuum case (i.e. J 1 = q = 0). 



The characteristic structure can be computed by considering the propagation of 
perturbations along a generic direction n l (belonging to an orthonormal tetrad 
{p,q,n}, where the index n stands for the longitudinal component). We consider 
an arbitrary solution (labeled generically by u) and formulate the eigenvalue problem 
for the perturbed fields (labeled by [u]). By defining /3 n = f3 l ni and computing the 
(Lagrangian) velocities v = v + j3 n , this eigenvalue problem can be written as 



$[&] = - e ijk nj [aB k ] + n'[a*] , (5) 

v[B l ] = e ijk nj [aE k ] + n^acf)] , (6) 

m = a[E n ], (7) 

m =a[B n ]. (8) 



It is now straightforward to obtain the following list of eigenvectors: 

• constraint modes: [<fi] ± [B n ] and [9] ± [E n ] , which involve the divergence-cleaning 
fields {<fi, ^} as well as the divergences of E and B. They propagate with light 
speed v = —(3 n ± a. 

• transversal modes: [E l — E n n l ] =p [e 8J rijBk], which can be written in tetrad 
components as [E p ] ± [B q ] and [E q ] =p [B p ]. They correspond to the EM waves 
and also propagate with light speed v — —f3 n ± a. 

Notice that due to the linearity of the Maxwell equations in vacuum, the eigenvectors 
depend only on the perturbed fields. Therefore, at this level, our discussion is 
completely generic and independent of the background solution u. 

Let us now consider a more realistic approach to our physical problem, -a spinning 
black hole interacting with a magnetic dominated, low density plasma, surrounding it-. 
The effect of this plasma is to provide both a charge q and current distribution J 1 which 
profoundly affect the behavior of the electromagnetic fields. As described in [9j 1], and 
summarized next, their role is accounted by explicitly defining the component of J 1 
transversal to B l and introducing a further constraint to account for the component 
of ,P parallel to B l , while q is obtained through Coulomb's law. In the following 
discussion, we will ignore the contribution of the divergence-cleaning terms since we 
are interested in describing the physical set-up. In our problem of interest, the plasma 
inertia is negligible, consequently the Lorentz force vanishes, defining the so-called 
force-free condition [T] , 

qE l + e^JjBu = . (9) 

This condition determines the current required to close Maxwell equations. The 
scalar and vector products of the force-free condition (j9|) with the magnetic field lead 
respectively to E l B t = and J 1 = qe ljk E 3 B k / B 2 + (J l B l )B l / ' B 2 . The component of 
the current parallel to the magnetic field can be obtained by imposing in the Maxwell 
equations (dt — Cp)(EiB t ) — 0, which is a natural consequence of the condition 
EiB 1 = 0. Additionally, in this electrodynamic limit the charge density is defined 
through the Maxwell constraint q — V^i? 4 , so that the current only depends on the 
electromagnetic fields, 

J 1 = (VmE^e^EjBk/B 2 + B l e^ km (B 3 V k B rn - E^ k E m )/B 2 . (10) 
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The force- free evolution system is obtained from the original Maxwell system (HE]) , by 
substituting current (flU)) in equation ([1]). Additionally, the field ^ and its associated 
evolution equation ([2J are eliminated since they become trivial with the definition of 
q = V ' iE l . The eigenvalue problem for the force- free approximation can be written as 

v[E l ] -- 

vm - 

v{4>] 

where S z 



- ae mk [B k ] + a {S l / B 2 )[E n ] + a (B l / , B 2 )j nk {B j [B k ] - E 3 [E k ]) , 
= ae mk [E k ]+an l [</>], 

= <*[B»]; (11) 

= e l i k EjB k is the Poynting vector. The diagonalization of the system is 
considerably more involved now due to the non-linearities introduced by the current 
pop . A convenient rearrangement of the fields makes this task easier, by changing to 
a basis containing the following combinations 





[C] 


= B q [Bp] — B p [B q ] 4 


E P [E q ] 


- E q [E p ] , 




[El + Bl 


= E p [E p ] + E q [E q } + 


B p [Bp] 


+ B q [B q ], 




[E t B l ] 


— E n [B n ] + B n [E n ] - 


f [EtB 1 ] 






{EtB*} 


= E p [Bp] + Bp[E p ] 4 


E q [B q ] 


+ B q [E q ]. 


The characteristic problem in this basis is given by 






( 


S n /B 2 


B n {B 2 - E 2 ) / B 2 
S n /B 2 





E n -1 






B n /B 2 















-E n _ 


1 





v U = a 




-E n B n /B 2 


-S n E n /B 2 1 





















1 












1 







\ 








E n 






(12) 



U 



(13) 



(14) 



with 

U=[C E n S n E 2 t +B 2 t 4> n B n EiB l ] T . 

The characteristic matrix can be diagonalized to obtain the following list of 
eigenvectors 

• constraint modes: [4>] ± [B n ], which propagate with light speed v = —j3 n ± a. 

• standing mode: [EiB 1 ] — E n [B n ], which propagates with speed v — — /3™, and 
contains information on [E^B ]. 

• "Poynting" modes: -B n A ± [C] + (S 2 A+ \~ -S n A ± )[ J B n ]- J B„ B n [B^+B^B 1 }, 
which propagate with speed v = —f3 n + qA ± , where we have defined A = 



±^VB 2 -E 2 . 



B 2 B 2 

• transversal radiative modes: B 2 (X + A" - l)E n [E n ]+2B n E n [E l B' l }+ B 2 (X + A~ + 
1 T 2S n ){ [E 2 + B 2 } ± [S n ] )T2E n B n [C] ± 2 E 2 n B n [0] , which are a generalization of 
the standard MHD Alfven modes, , and propagate with light speed v = — /3 n ± a. 

The force-free evolution system is strongly hyperbolic since there is a complete 
basis of eigenvectors for each direction ri 1 , ensuring the existence of a unique stable 
solution within the domain of dependence of the initial hypersurface in a boundary-free 
case. If the system is to be employed within a finite domain, boundary conditions have 
to be imposed on (some of) the fields, which may affect the stability of the solution. 
To ensure the well posedness of the resulting problem care must be taken to define 
these fields consistently One way to do so for symmetric hyperbolic systems, relies on 
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adopting maximally dissipative boundary conditions of the type U~ = R U + , where 
U T are the ingoing/outgoing eigenvectors of the evolution system, with \R\ < 1 [IT] . 
In a problem with constraints -as is the case here- further conditions must be derived 
to ensure the constraints are preserved. In the case of the force-free system as we have 
shown previously, the characteristic structure is rather complicated and depends on the 
background solution which complicates the analysis [J We defer to a future work the 
construction of constraint-preserving boundary conditions for the force-free system of 
equations. Here, since our main interest is to model different behaviours of the resistive 
load far from the black hole, we will consider the eigenvectors of the electrovacuum 
Maxwell system as a first step towards consistent and stable boundary conditions for 
the force-free system. Notice that this is short of defining conditions consistent with 
a force-free regime, however at far distances from the black hole, the plasma density 
decreases reducing to the vacuum case. Thus, depending on the physical scenario 
considered this choice will not be a physical limitation. Nevertheless, regardless of 
this issue, our approach will allow us to explore the robustness of the BZ mechanism 
by considereding different conditions and assesssing the stability and collimation of 
the resulting (if any) Poynting flux. The boundary conditions considered here take 
the form 

- E n n 1 } - [e ijk njB k }) = R ([E l ~ E n n 1 } + [e ijk nj B k ]) , (15) 

(M -[£»]) = R ([</>] + [B n ]) , [E n ]=0. (16) 

Notice that there is no rigorous proof that the resulting system will be stable with 
\R\ < 1 except in the electrovacuum case. Numerical experiments have shown that it 
is unstable if R < but stable otherwise. 

3. Results 

We adopt the formulation of the force-free approximation and the numerical 
implementation already described in [6j [7]. We will consider different boundary 
conditions corresponding to R = {0, 1/2, 1} in order to model the role of the resistive 
load far away from the source. Notice that we do not claim any of these values to 
be a true representative of the load -which is unknown in any case-, rather these 
values allow us to explore the resulting behaviour under profoundly different physical 
conditions and examine, in particular, the resulting Poynting flux behaviour. Indeed 
the resulting conditions contemplate that in absence of divergence errors, [B n ] = 0, 
and for: 

• R = 1: the transversal components of the electric field are free (i.e., determined 
only by the interior solution) while the transversal components of the magnetic 
field are set to 0. 

• R = 0: the transversal components of the magnetic field are equal of the 
transversal components of the electric field (and determined by the interior 
solution) . 

Any other value in R £ (0, 1) just spans the parameter space between these two limits. 

To examine the behaviour of the solution under these options we adopt a domain 
with (x,y) £ [—12, 12]M and z — [—16, 16]M covered with a uniform grid with grid- 
spacing given by A = 0.01M. We studied the three different cases evolving them 

| See e.g. [12| for a related example in the MHD case where constraint preserving boundary conditions 
arc defined. 
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Figure 1. Electromagnetic energy density flux for _R = (left), _R = 1/2 (middle) 
and R = 1 (right) at t = 320M. 



until a quasi-stationary regime is reached, comparing the luminosity obtained for 
each case far away from the black hole. This luminosity is computed as the integral 
of the electromagnetic energy flux density (which can be seen for the three cases 
studied in fig. H) ona solid angle of 15 degrees along the z-axis (further details on this 
computation are described in [7]) and shown in figure [5] (left plot), after a convenient 
renormalization with the asymptotic value of the R = 1 case. After a transitory initial 
behaviour, the obtained luminosities for the different values of the reflection coefficient 
R are qualitatively the same in structure and the emitted power within similar ranges. 
It is interesting to notice the energy density flux reaches a rather constant asymptotic 
configuration in the R — case, while relatively small amplitude cylindrical oscillations 
are observed for < R < 1. Nevertheless, the jet structure and average luminosities 
obtained are comparable in all cases. Thus, within typical astrophysical uncertainties 
it is clear the BZ mechanism is not affected when considering significantly different 
boundary conditions. 

To ensure the boundary location does not affect the above observations, we have 
studied the jet's luminosity for two different boundary placements along the jet's 
direction; the original domain (small) with z = [—16, 16]M and another one (large) 
with z — [— 32,32]M. Figure [2] (right plot) illustrates the obtained values for the 
case with R = 1 normalized to the asymptotic value of the small domain simulation, 
showing again that at late times the values obtained are quite similar. 

4. Final Comments 

This work illustrates the robustness of the BZ mechanism to generate a collimated 
Poynting flux of energy which is largely independent of the boundary conditions 
adopted. Consequently, while the load's characteristics are essentially unknown, our 
studies indicate the jet resulting from the plasma's ability to extract energy from the 
black hole's vicinity is robust. 
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Figure 2. (Left) Luminosity for the three cases considered, normalized to the 
asymptotic value of the R = 1 case. (Right) Luminosity for the two boundary 
locations considered in the case R = 1, normalized to the asymptotic value of the 
small domain simulation. 
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